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PREFACE 



The purpose of this research was to explore the relation 
between the horizontal wind field, the isobaric surface con- 
figuration and the vertical motion of the isobaric surface. 

This work was conducted at the U. 3. Naval Postgraduate 
School, Monterey, California, during the period January to 
May 195^. 

I wish to express my appreciation to Professor F. L. 
Martin for his advice and guidance throughout this study. 
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CHAPTER I 



INTRODUCTION 



It has long been recognized by meteorologists that there 
is an interdependence between the pressure distribution and 
wind flow on a level surface. In the steady state the pres- 
sure, Coriolis and centripetal accelerations all sum vector- 
ially to zero. The steady state is a condition which may be 
frequently approximated in nature, but it is probably seldom 
realized. Thus, there is a process of continual adjustment 
going on in the atmosphere. The wind accelerates in response 
to changes of pressure gradient, and pressure gradient changes 
in response to non-gradient winds. 

Rossby £ 8j and Cahn ^ 3 ] studied the mutual adjustment of 
pressure and velocity distributions in certain simple current 
systems. These currents were narrow but Infinite in length, 

.-ifc 

having a current profile similar to that which has now beoome 
associated with the Jet stream. 

Unbalanced forces appear to occur in the atmosphere in 
another important fashion, which has been emphasized recently 
by 0‘ Connor [ 6 ] because of its importance in prognostic anal- 
ysis. In this situation strong winds are found either upstream 
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or downstream of light winds. This condition is observed 
daily between the isotach maxima and minima occurring along 
the Jet stream on the upper-level facsimile charts prepared by 
the WBAN Analysis Center. Under these clroumstances it is easy 
to see that even though all forces may be balanced at a given 
instant, inertia may carry air into a region with speeds which 
are not compatible with the pressure distribution existing 
there. What occurs then is not well understood, but there 
must be some mutual adjustment between the wind velocity and 
the pressure distribution. 

Brunt [ 2 ] discusses an investigation of this problem by 
Brunt and Douglas. This treatment of the problem makes a num- 
ber of simplifications and does not show the effect of the 
curvature of the air motion. Fetterssen [?] discusses this 
problem in a slightly different manner, but makes the same 
approximations as Brunt and Douglas, The method of Petterssen 
is followed and extended in this study, but the equations of 

A 

motion will be expressed in terms of the topography of an iso- 
baric surface instead of the pressure distribution in a level 
surface . 

All equations and the discussion which follow refer spec- 
ifically to the Northern Hemisphere. 



2 



CHAPTER II 



DERIVATION OF THE ISALLOHYPTIC GRADIENT EQUATION 



The equations of frictionless horizontal flow on an ieo- 
barlc eurfaoe referred to a natural coordinate system are 



(f+KV)V -- 



and 



jX - 

dt 



a 






( 2 . 1 ) 



These may be written in vector form as follows: 

y /KX% 2 '(f-th V) Vt - 07 



( 2 . 2 ) 



Differentiating partially with respect to time, we obtain 

% + 1 /ff ) + V(K ^ + Vg) 



-nJL(dy) _ dV ^£7 
1 ^ J 7t ~dt 



Since 



M -ms 



and 



rt[(f+9KV+% )$ - VW^ KV ) +[ 

,/^ _ 2 -/d. 

V dt 



+ 07 







Taking a cross product of both sides with IX gives 



V Ttj ( 2 . 3 ) 
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+ 
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This is a form of the Brunt Douglas equation. The left 



aide, exceot for the factor c\ , is the negative of the isallo- 

0 

hyotic gradient. No terms have been ignored, and the only 
assumption thus far is that the motion is horizontal and fric- 
tionless . 

In this form (2.3) is difficult to evaluate, since every 
term on the right-hand side involves a time derivative. Other 
expressions may be obtained which may be substituted for some 
of these. One of these is obvious from the tangential component 
of (2.1). Thus may be substituted for UY in the first 

- ' '/T 

term on the right of (2.3). 

This type of direct substitution is not valid in the term 
A/ 

however, since the coordinate system used is not fixed 
(see Figure 1). In this case 



// 



77 






'? 



J . 



5 - SM 4 



and 



But 



Zt Wf / 



'n / i)P- -7 > 7 - o o \ 

_ -t >5 y 



’t \ 



GO 



z>x 



= . 






and 






700 {' 



Therefore, -j A/Kj 

f/tVir ' ~ 



•y- 



si? v -i-SO r ^V70 (JJ 

£ 7 ./ v / 4 r ' 



where (/'•' is the angle between the / axis of a fixed Zj coordinate 
system and the unit vector f which is directed along the 
positive 5 axis of the natural coordinate system. Throughout 
this work the 7 and “ axis of the fixed coordinate axis will 
be directed east and north, respectively. Thus, the angle V 
is the angle between the wind direction and the east. 
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Figure 1. 

Coordinate Systems Used in Derivation 



Differentiating 

jV 



av 



partially with respect to time gives 

29. - i>r 1 



or 



- 1 * *' V iU! v if + a ^ Ttj 



This may "be ‘written s _ . ■ ,, 

ty.\ - a.3- fir +«%£ &Z 

-ftr/T )- joi'-ot / * d On or . 

The expression on the right may then be substituted for 
) * n (2.3) and with the direct substitution for Sljf 

in the first term on the right side results in 

')2 ~r/ / 1 "5 u; i <3? . /W/r . , . ./ , . A -j\ 1 



cj V, $ -- - Jv[-(r+'H - Vg(f-n A V) + / V# 

+ j h ^ Vi 1 tvjy? ^ ii s ( : rt -] 



(2.4) 



Since the normal component of (2.1) may be written 

( / f HV )V i 'if?- 
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it la evident that the terra inside braces in (2.4) is zero. 
Also, the left-hand side of (2.4) may be written in component 
form as follows: 



The tangential components of the left and right sides of 
(2.4) are thus identical, and the right side adds nothing which 
will assist in evaluating the lsallohyptlc gradient. While 
there may be some method of evaluating the tangential compon- 
ent of (2.3) without a direct knowledge of the time rate of 
change of the acceleration of the wind, it is not obvious. In 
view of the fact that acceleration is usually ignored in atmos- 
pheric motion, it might appear that the term $■( j is negli- 
gible. If this is ignored, the tangential component of (2.3) 



is 









- V. 

> 



(2.5) 



A qualitative Interpretation of this equation with reference 
to the upper 1 air chart (Chapter III) is not too difficult, 



and an empirical check of it is discussed in Chapter IV. 
The normal component of (2.4) is 







( 2 . 6 ) 



Except for the term involving , all terms on the right- 

l' 

hand side of this equation can be evaluated from a synoptic 
chart. In order to explore the possibility of substituting 
another expression for , Blaton's equation relating tra- 
jectory curvature to streamline curvature is recalled. 



6 



It may be written h r * ■> 



/ ■ O' 



From this, taking the derivative with respect to time, 



v ; 

ot V 



y. 



and 






-/ - 



O' ( 



V 

J? 



(2.7) 



Now, 



3V 5 y’ _ v , 

V/ f " }h \ 

c - c, i ^ 



'IT i 



Carrying out the Indicated operation gives 



V< 









since 



r rjf 



This may be written 






"'Aj _ J 1 - ( -sry 

“Tf“ - 'c>5» ?? / •?/ ^ 

With this substitution, and after some rearrangement (2.7) 



becomes 



,V 






But 



and 



\ i | / .>■ i *> y , _ Sif— • — ! 

Sf u / + ' it < 



"n i/ 



1/ _ ’jj/ 

1)t Jt 



d '/ .. j * r 



\ A . A 

* v \S 



Therefore , v _ 

D\l / i / V , i n^i- if / \/ if - <=* A 

- Vi (tr / U v' * *r T$6t zs + 1 ^ Yr . 



1/ 



(2.8) 



Substituting this back in (2.6) yields 



* ssl 

,.-d «■/) it -j»f 



Collecting terms gives 







(2.9) 



The tangential component of the isallohyptic gradient may be 
taken in its complete form from (2.3). It is 



These two equations, (2.9) and (2.10), are the equations 
for the isallohyptic gradient normal to and along the wind 
vector, respectively. Equation (2.9) has been put in this form 
to make numerical evaluation of the equation from the weather 
map easier. There is some question whether this form is better 
than (2.6). It is to be noted that as written the units of 




( 2 . 10 ) 



Equations (2.9) and (2.10) are m sec 



8 



CHAPTER III 



A QUALITATIVE INTERPRETATION OF THE 
ISALLOHYPTIC GRADIENT EQUATIONS 



Although (2.9) represents the form In which the normal 
component is evaluated empirically (Chapter IV) , it is easier, 
perhaps, to gain an insight into the meaning of the equation 
from (2.6). 

^ / l) 

The term , if positive at a given point on an 

lsobarlc surface, Indicates that, relatively, the isobaric 
surface is rising more rapidly to the left of the wind direc- 
tion than to the right. Since the isobaric surface is lower 
to the left of the wind vector in virtually all instances 
above the friction layer, a positive value of this term indi- 
cates a relaxing lsohyptlc gradient. Similarly, a negative 
value of this term can be associated with a tightening iso- 
hyptic gradient. 

The tangential term of (2.10) is subject to a 

similar interpretation. If it is positive, the isobaric sur- 
face is rising more rapidly (or falling less rapidly) ahead 
of the wind vector than behind it. Conversely, if the term 
is negative, the isobaric surface is falling ahead of the wind 
vector and rising behind the wind vector relative to the point 
where the term is being evaluated. 

In order to evaluate (2.10), it is necessary to evaluate 
separately the terms (ffh and . The second term, 

t0 evaluate, since it is the local rate 
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of ohange of acceleration. It la quite small, probably never 
-7 

exceeding 10 m sec . On the other hand, the first term 
(-fiKV)V^f is also quite small. A method for evaluating the 
latter term numerically is discussed in Chapter IV; however, 
the sign of this term may be determined by inspection in par- 
ticular cases. 

The factors V and are virtually always positive 
in the Northern Hemisphere. Thus, the sign of the first term 
on the right-hand side of (2.10) is determined completely by 
the sign of the factor -~f . In a system of progressive open 
troughs and ridges such as that frequently found on upper- 
level charts, is positive at the ridge line and negative 
at the trough line. This indicates that on the basis of this 
term alone, the tangential component of the isallohyptic gradi- 
ent points in the direction of the wind at the trough line 
and opposite to the wind at the ridge line, provided the trough 



and ridge are progressive. 

■w 

r 

In order for the normal component to be positive, it is 
necessary that the sum of the three terms on the right-hand 
side of (2.6) be greater than zero. None of these terms are 
inherently either positive or negative, and each must be exa- 
mined carefully in a given situation to determine its sign. 
These terms will be taken up singly in the following paragraphs 
to try to determine approximate magnitudes and sign. 

Let /) : 1/-+^^ ^ be an abbreviated notation for 

the first term of (2.6). This term consists of three factors, 
of which is the only one which is always positive. The 
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factor is positive if the wind is blowing across isohypses 
toward higher values. The factor -f is always positive in 
the Northern Hemisphere. The wind velocity / is always 
positive. The curvature H way be positive (oyclonlo curvature) 
or negative (anticyclonic curvature). The term may be 
positive (counterclockwise turning) or negative (clockwise 
turning). The Coriolis parameter -f is of the same order of 
magnitude as %KV , though in most cases it is numerically 
greater. The term is usually not of the same order of 

magnitude as f , although there are occasional exceptions. 

From this one can ccnolude that the factor (f 
will usually be positive. An exception occurs when great 
anticyclonic curvature is associated with strong winds and 
the wind is locally turning clockwise with time. No such case 
has been found in the numerical evaluation of this factor. 

Thus, the sign of /\ is usually determined by the cross- 
contour flow of the wind. This term may vary in magnitude 

•it 

— 6 

from 0 to 10 m sec~' J . 

It is interesting to see what the contribution of this 
term is under various circumstances. Cross-contour flow to- 
ward higher values of Z gives a positive contribution to 

and indicates a weakening isohyptic gradient. Cross- 
contour flow toward low values of 2 indicates an Increasing 
Isohyptic gradient. The reverse of this would occur if 
(fi ?h were negative, which appears unlikely. 
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The second term, B~ , involves the factor 

. The variables appearing in this factor have Just 
been discussed. Except for the case of abnormal anticyclonic 
flow (V)>%Vgj , this factor is positive. The velocity,!/ , 
is positive always. Therefore, the sign of this entire term 
is the same as the sign of the factor -!-=• in nearly all cases. 

The sign of is positive if the wind increases along 
the streamline (i.e. weak winds upstream with strong winds 
downstream) and negative if the wind is decreasing along the 
streamline. The magnitude of the terra is a function of velo- 
city and rate of change of velocity along the trajectory in 
an obvious manner. It also is a function of curvature — 
strong cyclonic curvature being associated with greater magni- 
tudes, and strong anticyclonic curvature being associated with 
lowe r magni tude s . 

The term 8 alone then gives a positive contribution in 
almost every case when weak winds are moving into an area of 
strong winds.’ This is the entrance region of Scherhag C9 J • 
When such an entrance region is associated with cyclonic cur- 
vature, the contribution of B i s positive and may be relatively 
large. An entrance region associated with anticyclonic curva- 
ture will almost always give a small positive contribution. 

This difference between the cyclonic and anticyclonic cases 
i s again due to the factor of $ , discussed above. 

When strong winds are moving into a region of weak winds, 
the terra Q is nearly always negative. This is the delta 
region of Sherhag. A delta region associated with cyclonic 
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curvature results in relatively large negative contributions 
to . The contribution of the delta region associated 

with anticyclonic curvature is negative, but relatively srrall. 

In magnitude this term also may vary from 0 to 10“ m sec - ^, 
though it is doubtful that it approaches 10”^ very frequently. 

The third term on the right-hand side of (2.6) is C = • 

Thus, if the curvature is increasing locally with time, term C 
gives a negative contribution. This would occur in a region 
into which a trough is moving. Conversely, a ridge moving 
into a region would give a positive contribution to . 

The term C varies in magnitude from 0 to 10 - ^ m sec"*-^, 
and appears (See Chapter V) to be generally about one to two 
orders of magnitude smaller than term /] . The range of mag- 
nitudes given for all three terms must be considered approxi- 
mate only. 

Several model-t.voe examples are now presented in which, 
except as rioted, only the contributions due to terms // and B 
are considered. These examples are illustrated in Figure 2. 

The first example depicted by Figure 2(a) is a situation 
in which both terms A and B are positive. It should be 
noted that the winds, while flowing across contours toward 
higher values, are stronger downstream than upstream. There- 
fore, the lsallohyptic gradient associated with this must be 
negative (i.e. ^ is positive), and the isohyptic gra- 

dient is relaxing. 
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The curvature of the trajectories does not affect the 
sign of the isallohyptic gradient in this example unless the 
trajectory curvature is negative (anticyclonic) and exceeds 
the critical curvature for gradient flow. This case is dis- 
cussed below in connection with Figure 2(c). Trajectory cur- 
vature does change the magnitude of the isallohyptic gradient 
however, if all other things are held constant — cyclonic flow 
resulting in isallohyptic gradients of greater magnitude than 
anticyclonic flow. 

Figure 2(b) depicts a situation in which terms A and 3 
are both negative. The chief features to be noted are the 
flow of air across the contours toward the low and the decrease 
of wind speeds downstream. The isallohyptic gradient asso- 
ciated with this is positive and indicates that the isohyptic 
gradient is increasing with time. The influence of curvature 
is precisely the same in this case as for the first case de- 
scribed aboye . 

Figure 2(c) illustrates a condition often found at the 
crest of a ridge on upper-level charts. The curvature of 
the contours, given by R c , is negative, but of greater 
magnitude than the maximum curvature, given by Rj , for 
gradient flow with the existing contour gradient. For this 
reason, air which is moving over the crest of the ridge will 
be unable to make the sharp turn required to remain in the 
contour channel, and so will cross contours toward lower iso- 
baric heights downstream of the ridge. If the trajectory of 
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Figure 2. Schematic Illustration of Synoptic Situations for 
Which Sign of Isallohyptic Gradient May Be Determined by 

Inspection. 
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such an overshooting air particle is a trajectory of maximum 
curvature, then the factor -f-t'j h l//' is zero. The proof of 
this statement is given "below. 

Bjerknes fl j shows that for antlcyclonlc gradient flow 

/?„ = 

Now n - 

' St 

therefore, /?L - ~rr z - ^yr- 



Hence 

This condition probably is met occasionally in the situa- 
tion of Figure 2(c). Then K ~ /tw.jr and (2*6) becomes 



p-nist I - §5? 5s - V ^ 



( 3 . 1 ) 



In this case it is not immediately obvious that the 

2 \ i / ' 

term V ^ is small compared with the other terms of the 
equation, as has heretofore been assumed, nor are the two 



terms on the right side likely to be of like sign. Indeed, 
they are likely to give contributions of opposite sign. The 
factor will be positive, and 4J? will be negative, re- 
sulting in a negative contribution by the first term. The 
factor will determine the sign of the second term. It 
will generally be negative, resulting in a positive contri- 
bution by the second term. The sign of the lsallohyptic 
gradient, then, is determined by the sign of the term of 
greater magnitude. 



16 



This is a synoptic situation which meteorologists usually 
agree is accompanied by greater rises of height to the right 
of the air stream than to the left. These rises take the 
form of a propagation of the crest of the ridge downstream. 
Thus, it appears that the negative, or first term on the right 
of (3-D is the dominant one. In this example we may have one 
of the infrequent cases in which |Jr is of the same order of 
magnitude as f . 
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CHAPTER IV 



NUMERICAL EVALUATION OF THE 
ISALLOHYFTIC GRADIENT EQUATIONS 



Obtaining general solutions of (2.9) and (2.10) is not, 
in general, feasible. It is possible, however, to obtain 
numerical solutions in particular cases, subject to certain 
simplifying assumptions. 

Equation (2.9) may be written in terms of finite differ- 



ences along a trajectory, as follows: „ 

Mm - 



ac; 



v t . &% \ %\( 

At^\^r &ti{ ) ^ / \ 



h - % 

AV 



(4.1) 



The subscripts refer to points located along a series of 
trajectories as shown in Figure 3* 




Figure 3. Trajectories Showing the Grid Arrangement Required 
to Obtain a Numerical Solution 
of Equations (4.4) and (4.5). 
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Muoh of the success of evaluation depends on the accur- 
acy of the trajectories. Because of the fact that adequate 
upper-air charts are available only at 12-hour intervals, 
trajectories drawn by the method of successive approximation 
are of dubious value. Also, these require that the charts 
be available at both ends of the 12-hour time interval. This 
makes the result useless as a prognostic tool. It was de- 
cided to make use of constant absolute vorticity trajectories 
(CAVl^s) for the numerical work. 

The method of constructing CAVT’s from infle ctlon points 
is quite generally known. The assumptions Involved (non- 
divergence, constant lateral shear, frictionless flow, and no 
vertical motion) are quite restrictive. Equally restrictive 
for the purposes of this study is the requirement that the 
CAVT be constructed from an inflection point. Martin j^5 J 
in a study*. which was carried out concurrently with this study 
pointed out,how CAVT's might be constructed from any point. 

A modification of his method was used in this study and is 
described in Appendix 1. 

The use of a CAVT as the trajectory simplified the eval- 
uation of (4.1) in a number of ways. The arc length AS is 
equal to the product of V and At . The curvature , /f* , of 
the trajectory at point 2 (Figure 3 ) way be determined from 
the initial conditions and the latitude of point 2. Also, 
since the wind direction is known initially, it is a relatively 
simple matter to determine the turning of the wind at a given 
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point by measuring the angle between the initial wind direction 
and a tangent to the CAVT at that point. 

Equation (4.1) adapted to the use of 12-hour CAVT's be- 




(4.2) 



Equation (2.10) written in similar form for numerical 
evaluation is 



3&( n) - (HWK 



M, a 
Z /J 




(4.3) 



The second term on the right in (4.3) cannot be determined, 
but all other terms in (4.2) and (4.3) may be determined from 
a set of trajectories like those shown in Figure J>. Equations 
(4.2) and (4.3) with unit conversions, which permit data 
taken directly from an analyzed upper air chart to be used 



without change, become 

4- j«i h + 0 . 0 % w ) 



An 



+o, cn'Vi-v,)^ «'<it -HWi W- 



(4.4) 



and 



215 ^ 



(az) -- c.)% ( os-Ma td-)VA^ 



(4.5) 



The difference A A found on the left side of these 
equations is the 12-hour change of height of the isobaric 
surface. The units in which the various parameters are ex- 
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pressed, are as follows: \J (knots), K (reciprocal of degrees 

of latitude), £. (feet), ^ and (degrees), AV and AS 
(degrees of latitude). The lsallohyptic ascendant is expressed 
in terms of feet per degree of latitude per 12 hours. All 
computed and observed values are of the isallohyptic ascendant. 

A series of 44 computations were made using CAVT's on the 
500 mb. level charts for 0300? and 1200? , 12 March 1954 and 
03002 , 13 March 195^* The only criterion used in selecting 
the points from which the CAVT's were constructed was that 
initial conditions be as clear-cut as possible. Approximately 
one-third of the trajectories were started from troughs. One- 
third were from inflection points on the west side of a trough, 
and the other one-third were from inflection points east of 
a trough. The fourth term of the right side of (4.4) proved 
to be so small that it was ignored. 

Another series of 21 computations were made from CAVT's 
on 500 mb. level charts selected from the period 9 through 18 
April 195^* The criterion used in selecting the CAVT starting 
point was that it be on an isotach maximum. This made the 
second term on the right in (4.4) negative in every case. 

All trajectories of both series have starting points within 
the United States. 

The computed normal components of the isallohyptic ascen- 
dant are shown for the series of 44, with the observed values, 
in Table 1. The coefficient of correlation between computed 
and observed values is +0.65, Indicating that approximately 
43$ of the observed variation is accounted for by the computed 



21 



values. For 44 values the correlation coef f lclent , -f- 0 . 65 » 1 b 
highly significant at the 99$ level of belief, according to 
Hoel [ 4j. 

The regression equation (Figure 4 shows the scatter dia- 
gram and regression line of observed values if on computed 
values X ) is 

0^5 X . 

In this J and % are expressed in units of feet per degree 
latitude per 12 hours. This equation indicates that in this 
series, the computed values were, for the most part, of the 
right sign, but too great in magnitude. The dominant term 
in (4.4) was the cross-contour term in most of the computa- 
tions in this series. An attempt was made to reduce this 
term empirically by multiplying it by a constant factor. This 
did not improve the correlation, however. 

Table 1 also lists the computed and observed tangential 

.A 

components from this series of 44 trajectories. The correla- 

% 

tlon coefficient between these data is -f 0.59. indicating that 
approximately 35$ of the variance of the observed values is 
accounted for by the computed values. The correlation coef- 
ficient in this instance also is highly significant at the 
99$ level of belief. 

The regression equation of observed values on computed 
values is 

y = -5,3! +0,5/ X 

(See Figure 5 for the scatter diagram and plot of this equation.) 
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TABLE 1 



Components of Isallohyptic Ascendant Associated 
with CAVT's Constructed on a 5°0 Kb. Chart , 



Date 

and 

Time 


( 

Normal C 
Computed 


Isalloh 
feet per d 

omoonent 

Observed 


jptic Ascend 
eg lat per 

Tangential 

Computed 


ant 

12 hr a) 

Component 

Observed 


Remarks 




+ 53 


4 28 


0 

i 


- 18 






4 140 


-h 27 


- 46 


- 22 


From inflec- 




4 110 


4 21 


- 38 


~ 23 


tion point 




4 79 


0 


- 37 


- 26 


southeast to- 




4 3 


- 12 


- 8 


- 22 


ward trough 




4 120 


4 120 


0 


0 




03002 


- 9 


0 


0 


0 


From trough 


to 

1 


- 33 

! 


- 78 


0 


- 35 


northeast to- 


15002 


| "15 


- 26 


0 


- 4 


ward inflection 


12 Mar. 


t .32 


- 25 


0 

1 


- 4 


point 


1954 


4*24 

< 

— - ------ - . - - - - j 


- 29 


- 6 


0 






- 140 


0 


4 120 


4 33 






4 12 

' 


4 14 


- 20 


0 


From inflec- 




- 97 


1 - 10 

1 


4 110 


4 28 


tion point 




- 78 

1 


[ 

- 21 

1 


4 43 


+ 26 


1 

northeast to- 




— 44 


1 

- 12 


-f- 32 


4 9 


ward ridge 


15002 


4 170 


4 54 


- 51 


0 


' 


12 Mar. 


4 180 


4 91 


- 120 


- 37 


From inflec- 


to 03002 


4 110 


t 67 


- 75 


- 17 


tion point 


13 Mar. 


4 100 


4 38 


- 78 


- 33 


southeast to- 


1954 


f 30 


4 50 


- 14 


- 12 


ward trough 
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TABLE 1 
(Continued) 



Date 

and 

Time 


( 

Normal C 
Computed 


Isalloh 
feet per d 

omponent 

Observed 


yptic Ascend 
eg lat per 

Tangential 

Computed 


ant 

12 hrs) 

Component 

Observed 


Remarks 




+ 85 


-t 140 


- 40 


~ 46 






- 200 


- 91 


+ 140 


- 38 


From trough 




- 67 


- 45 


-t 15 


- 91 


northeast to- 


1500 ? 


+ 34 


- 5 


- 39 


- ^5 


ward inflection 


12 Mar. 


+ 27 


0 


- 36 


- 56 


point 


to 03002 


- 76 


- 62 


+ 38 


+ 4 5 




13 Mar . 


- 130 


- 43 


+ 97 


* 23 


From Inflec- 


1954 


- 65 


- 10 


+ 24 


+ 29 


tion point 




- 82 


0 


-h 3^ 


+■ 27 


northeast to- 




- 67 


-h 18 


-f- 41 


+ 31 


ward ridge 




+ 10 


+ 59 


+ 47 


+ 130 


From inflec- 


1 


- I 50 


+ 25 


+ 64 


+ 80 


tion point 




- 120 


0 


+ 57 


f 62 


southeast to- 




- 57 


- 6 


+ 16 


t 31 


ward trough 


0300 ? 


- 21 


+ 77 


- 83 


- 150 


From trough 


to 


-84 


- 37 


0 


- 27 


northeast to- 


1500 ? 


f 11 


- 25 


- 20 


- 42 


ward inflection 


13 Mar. 


-t 23 


- 5 


- 10 


- 18 


point 


195^ 


- 150 


- 72 


t 8 5 


f* 110 


From inflec- 




- 100 


- 16 


+ 38 


+ 42 


tion point 




- 85 


- 15 


f 39 


+ 31 


northeast to- 




- 53 


- 23 


f 44 


0 


ward ridge 




- 18 


__ 31 


0 


- 11 
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OBSERVED 




-200 -100 0 100 200 



COMPUTED 



Figure 4. Scatter Diagram and Regression Line of Computed 
against Observed Normal Isallohyptic Ascendant . 
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OBSERVED 




COMPUTED 



Figure 5* Scatter Diagram and Regression Line of Computed 
against Observed Tangential Isallobyptlc Ascendant. 
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TABLE 2 



Components of Isallohyptic Ascendant Associated 
with CAVT's Constructed from Isotach 
Maxima on a 500 Mb. Chart 



j 




Isallohyptio Ascendant 






(feet per deg lat per 


12 hr 8 ) 




Date 












and 


Normal Component 


Tangential Component 


Remarks* 


Time 

- 


Computed 


Observed 


Computed 


Observed 




195^ 

3/9-03 


- 40 


-5 


-h 4 


+ 17 


Infl. to tr. 


3/9-15 


- 170 


0 


+ 24 


0 


Near infl. to 












rdg. 


3/IO-03 


- 3 


+ 10 


0 


+ 29 


Rdg. to tr. 


3/10-15 


- 33 


- 40 


+ 34 


+ 17 


Infl. to tr. 


3/H-03 


- 210 


+ 38 


+ 160 


+ 43 


Tr. to. infl. 


3/H-15 


-f 28 


+ 14 


- 9 


+ 1^ 


Rdg. to. t<r. 


3/H-15 


- 280 


+ ^5 


+ 23 


+ 40 


Tr. arnd. top 


I 










of low 


3/12-03 


- 180 


- 100 


0 


+ 46 


Tr. arnd. top 


i 










of low 


3/12-15 


\ 

HV 

O 


0 


0 


+ 25 


Tr. arnd top 


1 










of low 


3/13-03 


- 160 


-hi 6 


+ 200 


+ 9 




4/10-15 


- 130 


-20 




- 33 


Tr. to. infl. 


4/11-03 


- 370 


0 


+ 230 


0 


Tr. to. infl. 


4/13-03 


■h 33 


- 8 1 


+ 13 


+ 46 


Infl. to. tr. 


4/14-15 


1 — 6 


+ 20 


+ 30 


- 17 


Tr. to. rdg. 


4/15-03 


f- 41 


- 23 


- 8 


- 14 


Infl. to. tr. 


V15-15 


+ 190 


- 17 


- 150 


— 54 


Near infl. to 










tr . 


4/16-03 


+ 150 


+ 150 


- 68 


- 19 


Infl. to tr. 
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TABLE 2 
( Continued) 



Isallobyptic Ascendant 
(feet per deg lat per 12 hrs) 



Date 












and 


Normal Component 


Tangential Component 


Remarks* 


Time 


Computed 


Observed 


Computed 


Observed 




1954 

4/16-15 


•f 180 


-56 


- 210 


-12 


Infl. to tr. 


4/17-03 


- 180 


f 140 


+ 110 


■f 12 


Near tr. to 










rdg. 


4/17-15 


t 95 


- 95 


-r 115 


+ 87 


Infl. acr. wk. 










rdg. and tr. 


4/18-03 


f 120 


- 40 


0 


0 


Wk. tr. acr. 



wk . rdg . and 
sir. tr. 



♦Abbreviations used in this column are: Infl. Inflection point, 

tr. trough, rdg. ridge, wk. weak, str. strong, arnd. around, 

acr. across, to. toward. 

% 
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As in the case of the normal component, this indicates that 
computed values were predominately of the right sign, but of 
too great magnitude. 

Of the set of 21 trajectories which were constructed from 
i so tach 'maxima , less than half resulted in computed normal 
components of the right sign. Also, as may be observed from 
Table 2, the magnitude of the computed values were almost 
Invariably larger than the magnitude of corresponding observed 
value s . 

The correlation coefficient between observed and computed 
values i s -0.09, which is not significant at any acceptable 
level of belief. 

The computed tangential components resulting from this 
set of 21 trajectories were somewhat better than the computed 
normal components. These (Table 2) were found to have a cor- 
relation coefficient of +0.41 vrith observed values. This is 
signif icant-tat the 93$ level of belief. Since it indicates 
that the computed values account for only 17$ of the variation 
occurring in the observed data, it cannot be considered a good 
forecast tool. In this case, as in the others, the magnitude 
of the computed values was too large. No regression equation 
was developed for either of the components computed from this 
series of trajectories. 

Because of the large contribution made by the cross- 
contour term in most of these computations, the CAVT's were 
examined to try to determine what might be done to the trajec- 
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tory to reduce this contribution. Usually, though not always, 
a reduction of the Initial wind speed would have resulted in 
a decrease in the contribution of the first term of (4.4), 
without changing appreciably the contribution due to the other 
terms. This might have improved the over-all results; however, 
more study is required before a definite statement can be made 
in this regard. 

In using the difference equations, (4.4) and (4.5) Instead 
of the differential equations, (2.9) and (2.10), it is assumed 
that the differences are small enough so that the difference 
equations closely approximate the Integrated forme of the dif- 
ferential equations. This assumption is almost surely not 
valid when the time interval is 12 hours. Even with a perfect 
12-hour trajectory, the results .would be somewhat uncertain 
because of this. 

Another source of error is found in the crudeness with 
which the basic data must be obtained. Curvature plays an im- 
portant part in this problem in several ways; yet it is diffi- 
cult to determine the curvature of the contours on an upper 
air chart with any degree of accuracy. The curvature of the 
air trajectory is even less attainable. Other parameters are 
equally difficult to evaluate. This is especially so in con- 
nection with short wave 8 for which the measurement of the 
Initial parameters becomes increasingly critical. This pre- 
sumably accounts for some of the lack of agreement of Table 2. 
Furthermore, very severe verification criteria were imposed: 
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The computed normal iaallohyptlc component was verified along 
a curve roughly orthogonal to the center point of the various 
CAVT’s. It is quite possible that only a slight displacement 
of the forecasted iaallohyptlc centers located by the computa- 
tional procedure would lead to the variation between the ob- 
served and computed iaallohyptlc gradients. The best test of 
the 12-hour forecast procedure would be to integrate the 
iaallohyptlc tendency equations to determine the complete fore- 
cast isallohyptic pattern, and then verify the location of rise 
and fall centers against the observed pattern. The forecasted 
isallohyptic gradients correlate sufficiently well (see Chapter 
V) with the observed gradients, in spite of the severe verifi- 
cation just mentioned. 

In certain of the cases presented in Table 2, the CAVT 
frequently was constructed from a minor trough or ridge or an 
inflection point between minor features. In such cases, the 
CAVT frequently crossed a synoptic feature such as a ridge (or 

•i. 

trough) , continuing right on part way across the next trough 
(or ridge) downstream. Trajectories of this type may be re- 
latively unrepresentative. The best test in suoh cases, also, 
would be that of verification of the forecast isallohyptic 
patterns . 
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CHAPTER V 



CONCLUSIONS 

In Chapter III it was shown that intensification of cer- 
tain synoptic models on an isobarlc surface can be determined 
qualitatively, provided the required parameters are determin- 
able. Of the parameters required, the most difficult to deter- 
mine are the cross-contour flow and the local time rate of 
turning of the wind. The cross-contour flow is by far the more 
important of these two parameters when using (2.6) or (2.9) 
with trajectories as drawn in this dissertation. 

O'Connor's hypothesis L 6j seems to indicate that the most 
important term in (2.6) is the second term on the right-hand 
side, since his discussion points out, primarily, those features 
associated with term g . This follows since O'Connor's point 
of view is that an isotach maximum or minimum will remain 
within its appropriate contour channel at each stage of its 
movement. Thus, no cross-contour flow occurs in this scheme. 
Furthermore, in many cases, O'Connor's method seems to give 
acceptable qualitative results. 

With the type of trajectories used in this work, the 
verification was also quite good, at least in the data of 
Table 1. As discussed in the next paragraph, the cross-contour 
term A is, in general, much larger than the longitudinal wind 
difference, term B , and the former gives an over-estimate 
(with correct sign) of the normal component of lsallohyptlc 



32 



gradient. Thus, the particular objective technique employed 
herein, while demonstrating the importance of term A , appears 
to over-estimate it and under-estimate term 6 . Both of these 
shortcomings might well be overcome by extracting data from a 
forecast map halfway through the forecast period. This would 
have required application of an iterative forecast procedure 
Involving time intervals of 1 - 2 hours at a step. However, 
the scope of the present work was to devise an objective fore- 
cast procedure using only the initial map. 

Some information concerning the relative importance of the 
various terms of (2.9) may be gained by further study of the 
more successful series of computations, the 44 which are sum- 
marized in Table 1. The average contribution of each of the 
terms (averaged without regard to sign) are: term A (cross- 

contour term) 73.4, term B (longitudinal wind difference 
term) 5-7* term C 2.2. The fourth term was ignored because 

it was small* compared with the other terms. The average ob- 

% 

served value is 34.8, and the average total computed value is 
78.1. Further, the coefficient of correlation between the 
values of the first term alone and the observed values is 
+0.55. This is significant at the 99$ level of belief. The 
coefficient of correlation between the second term values and 
the observed values is -0.03. This is not significant at any 
acceptable level of belief. Also not significant is the 
-0.04 correlation coefficient between the computed values of 
terms A and B . 
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Although terms d and C of (2.9) are small, there is 
some improvement in the correlation between observed and com- 
puted values when they are Included. The discussion at the 
beginning of this chapter seems to indicate that there must be 
a certain correspondence between the signs of the first and 
second terms. Thus, if the second term contribution can be 
used to Indicate the sign of the entire computed value, even 
though its value is only a small fraction of the total value, 
it should have the same sign as the contribution from the first 
or largest term. This is not borne out by the data of this 
series. Of the 44 pairs of values obtained by matching cor- 
responding terms A and B , there were 38 pairs having non- 
zero second term values. These 38 pairs consisted of 20 pairs 
having opposite sign and 18 pairs having the same sign. This 
lack of correspondence of signs is also indicated by the -0.04 
correlation coefficient between these values. 

In conclusion it appears that, while significant results 
have been obtained by using CAVT's drawn from selected points, 
some more reliable method of determining air trajectories is 
necessary to make this method a really useful prognostic tool. 
The evaluation of all terms of (2.9) and (2.10) is dependent 
on an accurate air trajectory, but it is in the evaluation of 
term A of (2.9) that trajectory errors are most keenly felt. 
Without actually constructing trajectories, however, some in- 
formation can be obtained by qualitative application of the 
principles which have been set forth. The examples discussed 
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in Chapter III represent synoptic patterns which are amenable 
to qualitative interpretation. 
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APPENDIX I 



A METHOD OF CONSTRUCTING CAVT ' S FROM ANY POINT 



After making the usual assumptions (i.e. non-divergence, 
constant lateral shear, horizontal frictionless motion) , the 
absolute vorticity equation may be written 



Vo 

~Z 



+ t = if 



r T 



( 1 ) 



The subscripts indicate initial conditions. 

Subject to the usual additional assumption that V is 
constant, /? may be determined at any latitude for a given set 

4 

of initial conditions. Thus 




or 



Substituting appropriate constants and conversion factors 

* 

to permit the use of V in knots and ff 0 and f( in degrees 
latitude changes this to 

D. V' 

‘ S/^(^n(p 0 -y<uni)i- X 

o 

Radii may be computed for the desired latitudes (either 
higher or lower than the initial latitude, depending on the 
initial wind direction) . Successive arcs of circles drawn 
end on end using appropriate radii result in the CAVT. A 
CAVT constructed by this means is not a curve of continuously 
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varying curvature. By increasing the number of arcs used in 
constructing the CAVT , it may, however, be made to approach 
as closely as is desired the "true” CAVT which has continuously 
varying curvature. Generally, the computation of radii at 
5 degree latitude intervals is adequate. 

Figure 6 illustrates the actual CAVT construction process. 
It is to be noted that the center of curvature of each succes- 
sive arc lies on a line drawn through the starting point of the 
aro and the center of curvature of the last arc. This insures 
continuity of the curve between arc segments. Table 3 lists 
the initial and computed data used in the CAVT construction 
shown in Figure 6. 



TABLE 3 

Initial and Computed Data Used in CAVT 
Construction Shown in Figure 6 



deg 


knots 


— SC 

deg 


R, 

deg lat 


i 


dfg 


+ t i 


R 

deg lat 


36.5 


85 


15 


6.0 


1*K2 


37.5 


13.8 


6.2 












^2.5 


11.7 


7.3 












^7-5 


9.7 


8.8 






i 


i 




52.5 


8.0 


10.6 
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Figure 6. Constructing a CAVT from a Point Other Than an 

Inflection Point. 
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